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ABSTRACT 


We investigate the possibility of perpendicular heating of electrons and the generation of 
" 90 °-electron conics" by particle diffusion in velocity space due to wave-particle interaction with 
intense auroral kilometric radiation. This interaction is made possible by the downward shift 
in the R-X cutoff below the electron cyclotron frequency that occurs in the presence of warm 
plasma. We simulate this condition and solve the diffusion equation using a finite difference 
algorithm. The results show strong perpendicular electron heating and indicate that the main 
characteristics of an electron conic distribution can be reproduced under these conditions. 
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I. INTRODUCTION 


The term "conic" was first applied to an electron distribution by Menietti and Burch 
[1985] and was used to refer to electron distributions that were enhanced along the loss cone. 
This type of conic distribution has been associated with upper hybrid resonance waves by Wong 
et al. [1988], suggesting that the electron distribution can be obliquely heated by this type of 
wave. Menietti et al. [1993] have observed another type of conic distribution involving 
perpendicularly heated electron distributions, which may be called 90°-conic distributions; these 
distributions were observed near auroral kilometric radiation (AKR) source regions. Subseque- 
ntly, Menietti et al. [1994] have hypothesized that the observed heating of electrons and the 
filling-in of observed loss cones is caused by particle diffusion in velocity space due to 
interaction with high intensity auroral kilometric radiation (AKR) (cf. their Figure 8). 

Auroral kilometric radiation is commonly observed at spectral densities of up to 10 -9 
(V/m) 2 /Hz [Gumett, 1974] and frequencies up to several hundred kilohertz and down to the R-X 
cutoff [Benson and Calvert, 1979; Gumett et al., 1983; Bahnsen, 1989]. Benson and Calvert 
[1979] find that the wave normal angle is probably close to 90°. In source regions, AKR is 
observed to exist at frequencies below the local electron cyclotron frequency [cf. Bahnsen et al. , 
1989; Menietti et al., 1993]. This occurrence is due to the presence of a relativistic plasma 
distribution, which can drive the R-X cutoff below the electron cyclotron frequency [Winglee, 
1985]. Ungstrup et al. [1990] have shown that the loss cone in the AKR source region is filled 
in and have hypothesized that this effect is due to fast diffusion brought about by a wave-particle 
interaction. In the present paper, we examine the hypothesis that fast electron diffusion due to 
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interaction of the electron distribution with AKR can cause the 90°-electron conic distributions 
observed by Dynamics Explorer, as proposed by Unstrup et al. [1990] and Menietti et al. 
[1994], 

Wu et al. [1981] approximate the diffusion coefficient for conditions applicable to 
fundamental R-X mode AKR. Using this approximation, we estimate the value of the diffusion 
coefficient and compute representative diffusion times for various distributions. To simulate the 
action of diffusion on a particle distribution, a finite difference algorithm is then used to solve 
the diffusion equation given by Wu et al. [1981] We shall show that perpendicular electron 
heating can be reproduced under conditions approximating those in the AKR source region. 

In Section II, we discuss some relevant observations. In Section III, we look at the 
phenomenon of diffusion and the effects of the form of the index of refraction and the resonant 
ellipse on that phenomenon. In Section IV, we simulate perpendicular diffusion due to AKR 
under plausible conditions and look at its effects on a model distribution function. In Section 
V, we draw conclusions concerning the likelihood of this mechanism. 
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II. OBSERVATIONS 


Plate 2 of Menietti et al. [1993] shows a spectrogram of an AKR event observed at an 
unusually high altitude in the nightside auroral region by DE 1 on 12 October 1981. Two 
spacecraft positions where intense AKR emission extends below the local electron cyclotron 
frequency are indicated on the spectrogram, one at 00:28:30 universal time (UT) and the other 
at 00:29:30 UT, which Menietti et al. [1993] (cf. , Bahnsen et al., 1987, 1989) conclude are 
source field lines. In this plate, a region is indicated between two dashed lines where strong 
perpendicular heating was observed. This region extends from poleward of the observed source 
regions to a position between the two source regions. 

Figure 1 of the present paper shows a contour plot of a perpendicularly heated electron 
distribution occurring adjacent to and poleward of the source region observed in the 
aforementioned event at 00:28:30 UT. This figure shows a distinct enhancement of the electron 
distribution perpendicular to the magnetic field extending to about 4.5 keV. The electron 
distribution in the parallel direction, which we assume is unaffected by particle diffusion, extends 
to about 0.6 keV. There is also a clear filling in of the loss cone. We believe that the observed 
asymmetry and filled-in loss cone are the result of a wave-particle interaction involving intense 
AKR. We shall attempt to reproduce these effects by modeling the radiation spectrum, the 
electron distribution, and the diffusion process. 
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III. THE DIFFUSION COEFFICIENT, INDEX OF REFRACTION, 
AND RESONANT ELLIPSE 


The diffusion coefficient for the case of R-X mode waves propagating nearly 
perpendicular to the magnetic field is given by Wu et al. [1981] as 

D i± (v A ,v.) = [ d 3 k 5(t) - Q c h - k |V| ] (1) 

where ^ k (t) is wave energy per mode given here as a function of time, k is the wavevector, u r 
is the real part of the wave frequency, 7 = l/v/(l-v x 2 /c 2 - v^/c 2 ), and fi c is the absolute 
value of the electron cyclotron frequency. For | k | v j |/a) r <? 1, the evolution of the electron 
distribution function is given by the following equation from Wu et al. [1981] 



z- 4- D(v ^ v i> ^ 


v i 9v i 


( 2 ) 


In this study, we do not consider the effect of the wave-particle interaction on the wave 
intensity, so that the wave energy per mode is considered to be constant. 

Setting the argument of the delta function in Equation 1 equal to zero for a given plasma 
frequency, cyclotron frequency, and wave normal angle, yields the resonant ellipse, which is the 
locus of points in velocity space that can interact with the wave at a given frequency. Particles 
not on a resonant ellipse do not interact with the wave. The parameters of the resonant ellipse 
are given by Melrose et al. [1982] as 
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where C, is the displacement of the center of the ellipse along the vj/c axis, a is the semimajor 
axis, and b is the semiminor axis. All parameters given here are written as a fraction of the 
vacuum speed of light. 

Benson and Calvert [1979] determined observationally that AKR has a low-frequency 
cutoff at the R-X cutoff frequency. However, as mentioned above, Viking and DE 1 
observations frequently show that the low-frequency cutoff of AKR can reach below the electron 
cyclotron frequency in the source region, indicating that the R-X cutoff has been lowered. The 
lowered R-X cutoff is an effect of the presence of relativistic electrons, as shown by Winglee 
[1985]. This effect is crucial for the existence of perpendicular electron heating by AKR. 

In Figure 2 we plot the cold plasma index of refraction, in the top panel, and the square 
of the semimajor axis of the resonant ellipse, in the bottom panel, as a function of frequency for 
the following values of cyclotron frequency, plasma frequency, and wave normal angle, 
respectively: f c = 60 kHz, f p = 6 kHz, 0 = 105°. The R-X cutoff frequency is f,^ = 60.6 
kHz; this frequency is the lower limit of the abscissa. It can be seen that the index of refraction 
starts at zero at the R-X cutoff and rises asymptotically toward one. The bottom panel shows 
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that there is a limited frequency range for which the square of the semimajor axis has positive 
values, starting at f = 60.8 kHz and continuing to higher frequencies. Resonant ellipses do not 
exist outside this frequency range for the R-X mode, at the specified plasma frequency, 
cyclotron frequency, and wave normal angle. Figure 3 shows several of the resonant ellipses 
contained in this region. Note that the size of the ellipses decreases to zero as the frequency 
approaches its limit at either end of this range. Note also that the area where particle diffusion 
is allowed is limited to the region in which the ellipses occur, to one side of the origin, as 
expected from Wu [1985]. 

Because none of the resonant ellipses contain the region where vj/c » 0, it is impossible 
to produce perpendicular heating. However, as shown in Figure 1 of Winglee [1985], the 
presence of warm plasma lowers the value of the R-X cutoff frequency. Equations 4 and 5 show 
that lowering the value of a> increases the size of the resonant ellipse. Thus it will prove useful 
to study the effect of lowered wave frequencies on the configuration of resonant ellipses. 

A comprehensive warm plasma calculation [cf. Winglee, 1985] is beyond the scope of 
this paper. Therefore, in order to simulate the relativistic plasma effect on the R-X cutoff and 
maintain the simplicity of the cold plasma approximation, we have simply evaluated each point 
of the cold plasma index of refraction for a new frequency shifted downward by an arbitrary 
factor A times (f^ - f c ), where f RX is determined by cold plasma theory. Figure 4 shows that 
the result of this operation is a reasonable approximation to the actual warm plasma index of 
refraction. In this figure we replot data from Figure lb of Winglee [1985] for a DGH 
distribution function with v Xh = 0.3c in terms of index of refraction as a function of wave 
frequency. We compare this result to the cold plasma index of refraction curve computed from 
the same plasma frequency, cyclotron frequency, and wave normal angle and new transformed 
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frequency until n = 0 at f = f^, evaluated from Equation A/2 of Winglee [1985]. Clearly, the 
two curves are similar in shape, rising abruptly from n = 0 at the cutoff and then approaching 
n = 1 asymptotically. The only visible difference between the two curves is that the value of 
dn/dw near the cutoff is larger when calculated from the warm plasma distribution than for the 
frequency-shifted cold plasma distribution. This simple approximation is not rigorously correct, 
but it does preserve the important qualitative features of the true index of refraction and is 
analytically very simple. 

The effect of the downward frequency shift on the allowed resonant ellipses is shown in 
Figures 5 and 6. Figure 5 shows, in the top panel, the index of refraction, identical to that in 
Figure 3, shifted downward in frequency by 2^ - Q- In the bottom panel, the effect on the 
square of the semimajor axis is displayed. Two effects on the resonant ellipse are obvious: 
first, the range in which the square of the semimajor axis is positive is extended, beginning at 
the new R-X cutoff frequency and extending to about the same upper limit as shown in Figure 
2; second, the semimajor axis starts at a nonzero value at the R-X cutoff (fRx = 59.4 kHz) and 
generally has much larger values than in Figure 2. The semimajor axis is seen to be nonzero 
at f = f^, where n = 0. The resonant ellipse corresponding to f = f^ is centered on the 
origin, according to Equation 3. Figure 6 displays some of the ellipses in the frequency range 
where the square of the semimajor axis is positive, as shown in Figure 5. This simple exercise 
has shown that merely shifting the index of refraction to lower frequencies can cause the 
resonant ellipses to intersect the region surrounding the origin of velocity space, a necessary 
condition for a 90° -electron conic distribution to be produced. 

Figure 7 shows five cuts of the function D x i(Vx,V|) for a = 1 kHz, A = 2.0, and 6 
= 105° for different values of vj. It is seen that the values of D ± A are on the order of 10 17 
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cm 2 /s 3 . The sharp decrease in the value of D x x to zero occurs when there is no zero of the 
argument of the delta function in Equation 1, i.e., when 


2 v 2 

>w r -k|(co r )V| 

c 2 c z 


( 6 ) 


for all possible oo r . 

As shown by Figure 4, cold plasma theory yields an underestimate of the value of the 
value of dn/du, implying an overestimate of the group velocity. A lower group velocity should 
allow the wave more time to interact with the electrons near the source region and thus result 
in more efficient heating of the plasma. 
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IV. DIFFUSION MODEL 

We make two major assumptions in order to simplify the calculation of Equation 1. 
First, we assume that the AKR is emitted at a single wave normal angle. Here, f c = 60 kHz 
and f p = 6 kHz for all cases, similar to observed values from Menietti et al. [1993]. Energy 
density is given by the integral of the energy per mode, as 

W = ejEl 2 = jd 3 k % W 

where ^ is the energy per mode. To describe a wave impinging on the plasma from a discrete 
direction, we write the energy per mode as 

% = K(k) 5(6 - 6 0 ) 8(<(> - 4> 0 ) < 8 > 

k 2 sin0 

where K(k) is a function to be determined. The expression for the energy density becomes 

W = e 0 1 E | 2 = fdkK(k) 

J (9) 

= j du ^ K(k) . 


The energy density can also be written as 



By comparing the integrands of the last two equations, 


™ ^ 


( 10 ) 


( 11 ) 


By inserting this expression into Equation (1), a simple expression can be obtained for D Xi , 
evaluated at 6 = 6 0 , 4> = <t> 0 : 
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D i±( v x’ v i) = 4l ^ r { dc0 e ° d ^ ~~ 5 ^ r _ n c h - k i v |) • 

The value of d | E 1 2 /doo can be obtained from the power spectral density of AKR. A 
typical peak value of this quantity is 10~ 9 (V/m) 2 /Hz. Since the spectral density of AKR goes 
to zero at f R . x , reaches a peak, and then falls off rapidly with increasing frequency, we have 
approximated the spectral density by the form 

f ~ fRX exp{-[ f ~ f H] 2 } (13) 

a a 

where f > fRx and a is a free parameter. Figure 8 shows this function for A = 2 for five 
values of a. Figure 9 shows this function for A = 2 and a = 200, 500, and 1000 Hz. 

The assumption of a discrete wave normal angle for AKR emission makes the 
computation of the diffusion coefficient by way of Equation 1 straightforward. One effect of 
this assumption is that there is an abrupt break between zero and nonzero values of the diffusion 
coefficient corresponding to the distinction between points on and not on a resonant ellipse. The 
discontinuity is a nonphysical result of the assumption that the radiation occurs at a single wave 
normal angle. We have artificially removed the sharp edge by smoothing the results. 

We now show that the model assumed for the dispersion relation in the last section can 
be used to generate conic-like distributions for an electron population in the presence of intense 
AKR. We assume the electron population is physically adjacent to the AKR source region, as 
observed, and can be described by two components. The first component is a cold Maxwellian 
of energy 10 eV. The second component is a "ring distribution" with radius of v 0 and a loss 
cone of Ioss . The ring distribution is calculated from the formula 
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( 14 ) 


where a loss = 40° and v 0 = 0.075 c. A contour plot of this initial distribution is shown in 
Figure 9. This model does not include a relativistic population of electrons sufficient to lower 
the R-X cutoff significantly. According to Winglee [1985], such a distribution must have a 
thermal energy around 0.3c to be effective. The cool distribution shown in Figure 9 should be 
considered to be superimposed on the dominant relativistic plasma. 

With the stated conditions, we have solved Equation 2 with a finite difference algorithm 
for two different sets of parameters. The time increment is 10* 3 T diff , where r diff is given by r diff 
= v lh 2 /D' L lt D' x A is taken as the value at the origin of velocity space, and v th is the thermal 
velocity of the Maxwellian component of the unperturbed distribution function. For the first 
case, a = 300 Hz, A = 1.5, and 0 = 100°. The result of running the finite difference 
computation for three diffusion times are seen in Figure 10. In comparing the result with Figure 
1, we see that the Maxwellian portion of the model distribution is stretched similar to the data, 
although not quite as far. The loss cones are seen to fill in. In Figure 9, the sharp edge at vj/c 
= 0 and at v x /c ~ 0. 10 is attributable to a sharp edge in the diffusion coefficient as discussed 
above. 

Figure 11 shows a case with parameters a = 1 kHz, A = 3.0, and 0 = 105°. The 
running time in this case is also three diffusion times. The Maxwellian portion of the 
distribution function is stretched and the loss cones filled in, similarly to Figure 10. The 


diffusion "edge" is not observed. 
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Figure 1 differs from both Figures 10 and 11 in that both the lower and higher energies 
of the distribution function appear to be unaffected by diffusion. We shall discuss possible 
reasons for this discrepancy in the final section. 
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V. CONCLUSIONS 

This work has shown that electron diffusion due to particle interaction with strong AKR 
is an effective perpendicular heating mechanism for electrons. In particular, this mechanism is 
capable of generating "90°-electron conic" distributions on short time scales in regions adjacent 
to the AKR source region. The essential physics that enables this process to occur is shown in 
Figures 2 through 6. These figures show the effect of warm plasma in lowering the value of 
the R-X cutoff frequency, as discussed by Winglee [1985], enabling the resonant ellipses to 
encompass the origin of velocity space. Without the presence of warm plasma, the resonant 
ellipses would all lie to one side of the v x axis of velocity space, as shown by Wu [1985], and 
perpendicular diffusion due to AKR would have little effect on the electron distribution. 

The results show a good qualitative agreement between the data and the model. In 
particular, we note a perpendicular stretching of the distribution function combined with some 
filling of the loss cone. However, we also note some significant differences between the data 
and the simulation. (1) The data shown in Figure 1 indicate that heating occurs most efficiently 
at intermediate energies, leaving low and high energy components of the distribution nearly 
unchanged. We were unable to recover this effect. (2) Our computation gives a sharp transition 
between zero and nonzero regions for the diffusion coefficient in the v x - V| plane. This edge 
shows up in the simulation as a discontinuity in the evolved distribution visible in Figure 10. 

Clearly, the approach to heating by diffusion undertaken here is only adequate to give 
a qualitative result. The weaknesses inherent in this approach are (1) use of a Maxwellian in 
the initial distribution since the initial distribution shown by the data is clearly non-Maxwellian; 
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(2) the downward frequency shift used to simulate the effect of a relativistic particle distribution; 

(3) the use of a simple spectral density distribution as given in Equation 9 an Figure 9; and 
especially (4) the use of a single wave normal angle to describe the radiation distribution. 
Introducing a range of wave normal angles in Equation 1 would tend to soften the sharp 
transition edge of the diffusion coefficient between zero and nonzero values. 

Perpendicular heating of electrons due to AKR is inherently a warm plasma effect and 
is expected in regions adjacent to the AKR source region. The generation of 90°-electron conics 
near the AKR source, typically at altitudes of 3000 to 5000 km in the nightside auroral region 
can lead to adiabatic folding of the distribution and more typical, smaller pitch-angle, electron 
conics observed at higher altitudes by DE 1 [cf. Menietti et al., 1994, Figure 8]. 
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Figure Captions 

Figure 1. Contour plot of the electron distribution function on day 285 at 00:28:04 UT, 
poleward of an AKR source region. This contour shows perpendicular heating 
of the low temperature distribution out to about 4.5 keV. (From Menietti et al. 
[1993].) 

Figure 2. (Top) The cold plasma index of refraction plotted against frequency for 0 O = 
105° , f ce = 60 kHz, and f pe = 6 kHz. The lower frequency limit, about 60.59 
kHz, is the R-X cutoff. (Bottom) The corresponding values of the square of the 
semimajor axis of the resonant ellipse plotted against frequency. This figure 
shows that there is a limited range of frequencies in which the resonant ellipse 
exists. 

Figure 3. Resonant ellipses corresponding to the region of positive semimajor axis squared 
shown in the previous figure. Note that the range of vj/c is from -0.1 to -0.3. 
Thus, none of these ellipses encompass the origin. 

Figure 4. Comparison of a numerically calculated warm plasma index of refraction, from 
Winglee [1985], for a DGH distribution with a = 1, v-n, = 0.3c, fi c = 60 kHz, 
u) p = 27 kHz , and 0 = 80° with a cold plasma distribution function of the same 
fl c , dip, and 0 transposed downward in frequency so that the cutoffs coincide. 

Figure 5. (Top) Index of refraction, identical to that shown in Figure 2 but transposed 
downward in frequency so that the R-X cutoff is below the electron cyclotron 
frequency. (Bottom) Square of the semimajor of the resonant ellipse for the 
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Figure 6. 


Figure 7. 


Figure 8. 


Figure 9. 


Figure 10. 


Figure 11. 


index of refraction shown in the top panel. Compare with the bottom panel of 
Figure 2 to see how the transposition of the index of refraction downward in 
frequency changes the value of this quantity. 

Resonant ellipses corresponding to the region of positive semimajor axis shown 
in the previous figure. 

Values of D ii (v i ,V|) as calculated from Equation 1 for fixed values of vj. 
Plasma parameters are f p = 6 kHz, f c = 60 kHz and wave normal angle 105°. 
The distortion parameter is A = 2.0 and the spectral width a = 300 Hz. 

The function, given in Equation 14, used to model the power spectrum of AKR. 
The modeling is done for five values of a. 

The initial distribution function used for the simulation run. This initial 
distribution consists of a Maxwellian of thermal energy 10 eV added to a ring 
distribution of the form given in Equation 7 with v c = 0.075c and of loss = 40°. 
The result of evolution of the initial distribution function shown in the previous 
figure after evolution for 3r diff , where r dlff = 0.27 s. The maximum of the v x 
scale is approximately the same as that in Figure 1. In this simulation, a = 300 
Hz, A = 1.5, and 6 = 100°. 

Same as Figure 9 but for a — 1000 Hz, A = 3.0, and 6 = 105°. The evolution 
time is approximately 3r dirf , where r diff = 0.34 s. 
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Contour Plot of Distribution Function 





